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Abstract: In this contribution we summarize experimental information and theoretical results for 
the dissociation cross sections of charmonium by light hadrons. Theoretical predictions for these 
RHIC-relatcd processes differ by orders of magnitude over the physically relevant energy range. The 
results found by the author and collaborators using a constituent interchange model, which predicts 
cross sections in the mb region near threshold, are discussed in more detail. 



1. Introduction 



Many unusual subjects have been studied in the 
name of QCD. One of the more unusual, which 
has arisen in the field of heavy ion collisions, is 
the size of cross sections of charmonia on light 
hadrons. This has attracted attention because of 
the proposal by Matsui and Satz |Q that suppres- 
sion of J/tjj production could be used as a signa- 
ture for the formation of a quark-gluon plasma. 

This suggestion, like many signatures pro- 
posed for the quark-gluon plasma, is perhaps ex- 
cessively intuitive. The idea is that a QGP will 
screen the linear confining interaction between 
quarks, so that a cc pair produced within a QGP 
will be less likely to form a bound cc charmonium 
resonance, as in Fig.l, but instead will more likely 
separate to form open-charm mesons. 

Even if this simple picture of cc production 
in a QGP is qualitatively correct, it can only be 
confirmed easily if the competing direct charm 
production and scattering by the initial relativis- 
tic nucleons is understood || and if there is little 
subsequent dissociation of the charmonia by the 
many other "comoving" light hadrons produced 
in such a collision. To summarize the last point, 
if charmonium + light hadron "comover" disso- 




Figure 1: The two scenarios we wish to distinguish 
for the external evolution of a J/tp produced (pur- 
portedly with reduced probability) within a quark- 
gluon plasma: 1) weak, versus 2) strong J/ip absorp- 
tion by comoving light hadrons. 
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ciation cross sections are small (case 1, top of 
Fig.l) and the background of direct charm pro- 
duction from the initial nucleons is understood, 
one may have a useful signature of QGP forma- 
tion, but if the comover dissociation cross sec- 
tions are large (case 2, bottom of Fig.l) one must 
distinguish a QGP-reduced charmonium produc- 
tion amplitude from subsequent dissociative scat- 
tering, and the interpretation of this signal will 
therefore be ambiguous. 

Thus it is of great relevance to the interpre- 
tation of RHIC physics to establish the approxi- 
mate size of these low energy cc + light hadron 
cross sections. 



where p is the mean nucleon density and L is 
the estimated mean path length in the exper- 
imental nuclear system. A "naive" interpreta- 
tion of the J/ip production data from collisions 
of various nuclear species using this formula gives 
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6 mb at yfs « 10 GeV [§, with a nu- 
merically similar result for the ip'. Of course 
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2. Experiment, or what passes for it 

Unfortunately we have no charmonium beams or 
targets, so the experimental cross sections must 
be inferred indirectly and are poorly known. The 
earliest estimates of lower energy charmonium 
hadronic cross sections came from J/ip photo- 
production experiments in the mid 1970s, which 
were interpreted in terms of a J/ip + N cross 
section given additional theoretical assumptions. 
Early Fermilab and SLAC photoproduction ex- 
periments gave rough estimates of ~ 1 mb for 
cj/ip+Ni assuming vector dominance, for photon 
energies from E 1 « 13 to 200 GeV |, §. A 
subsequent SLAC photoproduction experiment 
in 1977 used the A dependence of J/ip absorption 
to infer a rather larger cross section of o,j/^+n = 
3.5(0.8) mb at E 1 « 17 GeV (*/s w 6 GeV) §. 
The vector dominance hypothesis may have lead 
to an underestimate of the cross section in the 
earlier references ||. 

In heavy ion collisions these cross sections 
may be estimated from the ratio of lepton pairs 
produced in the J / ip peak to "background" Drell- 
Yan pairs nearby in energy. Since the J/ip must 
reach the exterior of the nuclear target to decay 
into a sharp mass peak, this ratio gives us an 
estimate of the absorption cross section through 
a classical survival probability formula, 

cr(J/ip — ► u + fi~)/cr(Drell — Y an a + fjT) 
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Figure 2: A fit of (2.1) to experimental J/ip pro- 
duction versus path length jjj; the line corresponds 
to 6.2 mb. 

one may raise many questions about the validity 
of this simple estimate, such as the importance 
of shadowing in Drell-Yan, the neglect of J/tp 
scattering by other light hadrons formed in the 
collision (such as ir and p), and the assumption 
of a single, constant J/tp + N cross section in all 
circumstances. 

Recently, concerns have been expressed that 
the J j ip and ip' wavefunctions have not had suffi- 
cient time to form within the nucleus in these col- 
lisions, so experiment may instead be measuring 
the cross section for a small initial cc "premeson" 
on a nucleon. One can increase the time spent 
in the interior of the nuclear system by selecting 
small and even negative xf events, as has been 
done by E866 at Fermilab. As discussed by He, 
Hiifncr and Kopeliovich H, this leads one to in- 



fpr rr abs 



2.8(0.3) mb and af, s +N = 10.5(3.6) 



exp(-p a a j% +N L) 



(2.1) 



mb respectively, also at w 10 GeV. This is 
rather more satisfying to people who have an in- 
tuitive notion that the larger -0' should have a 
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larger reaction cross section. Actually the con- 
nection between cross section and physical extent 
is less direct (compare K N and KN) , and in any 
case the relative proximity of inelastic thresh- 
olds alone would suggest a larger ip' cross section. 
These experiments also indicate a preference for 
dissociation over elastic cross sections in this en- 
ergy region by roughly a factor of 30 ||| . 

3. Theory: Introduction and pQCD 

To quote B.Miiller, "...the state of the theory of 
interactions between J /ib and light hadrons is 
embarrassing []. Only three serious calculations 
exist (after more than 10 years of intense discus- 
sions about this issue!) and their results differ by 
at least two orders of magnitude in the relevant 
energy range []. There is a lot to do for those 
who would like to make a serious contribution to 
an important topic." 

The theoretical situation has improved in the 
subsequent year, at least in terms of the quantity 
of calculations if not in the understanding of the 
scattering mechanism. A list of cc + light hadron 
cross section calculations is given in Table I. 

The most cited work, albeit furthest in its 
predictions from a low-energy "theoretical mean" , 
is the pQCD calculation of Kharzeev and Satz 
]lC| . This reference is basically a restatement 
of the color-dipole scattering model developed 
in the late 1970s by Peskin and Bhanot |ll|] for 
scattering of light hadrons by Coulombic bound 
states of very massive quarks. According to Pe- 
skin, the criterion for validity of this approach is 
"...not met even for the bb system." so there 
may be large systematic errors at the cc mass 
scale. This approach certainly makes marginal 
approximations for charmonium, such as the use 
of Coulombic wavefunctions (which are far from 
accurate for cc) and the introduction of a QQ 
binding energy (which is hard to interpret for 
charmonium, and is taken to be 2 Mtj — Mj/^, by 
Kharzeev and Satz). These color-dipole scatter- 
ing formulas also implicitly assume that charmo- 
nia are small relative to the natural QCD length 
scale. Since potential models actually find rms cc 
separations of about 0.4 fm for the J/ip, 0.6 fm 
for the Xc states and 0.8 fm for the ib 1 (lD< this 
is also a dubious assumption. 



Method 


Init. State 


Ref. 


pQCD 


J/ip + N 


1 


10 




meson ex. 


J/lp + TT 
J lib + TT, P 
J/lp + TT,p 

J/ib + N 

J/ib + TT, K, p, N 


1 
1 

1 
1 
I 

k 


H 

1C 
17 
, 


1 

L9| 


constit. int. 


J /lb + TT 

J jib + tt, p ; ib' + tt,p 
J/i{j + tt,N ; V' + tt,N 


fcf 

i 
i 


21 
22 


2d 



Table 1: A summary of cc + light hadron cross 
section calculations. 

Although this approach has problems with 
justification for cc, the predictions are nonethe- 
less interesting as estimates of the scale of these 
cross sections assuming a color-dipole scattering 
mechanism, and the approach could presumably 
be extended to lower energies by generalizing the 
wavefunctions and interaction. The formula for 
the J jib + N cross section quoted by Kharzeev 
and Satz [|10| is 



a 



J/ip+N 



2.5 mb 1 



Ao 
A 



6.5 



(3.1) 



where A = (s - Mj,^ - M%)/2Mj/^, the con- 
stant Ao is defined to be " ~ Mn + eo" according 
to the text following Eq.(24) of jio) (we assume 
the equality), and the "binding energy" eo is set 
to 2M D — Mj/^p. The result is shown in Fig. 3, 
together with the single lower-energy SLAC ex- 
perimental point j^j. 

Evidently the Kharzeev-Satz cross section is 
smaller than this SLAC point (which was an in- 
ferred cross section and certainly needs confir- 
mation) by about two orders of magnitude, and 
falls precipitously as ^/s is decreased. This ap- 
proach actually has no direct information about 
physical thresholds, so Kharzeev et al. typically 
leave their curves "dangling" just below y/s = 5 
GeV. (See Fig.2 of @ for example.) If we actu- 
ally plot their formula (3.1) for J/ip + N at low 
energy (Fig. 3), we find the unphysical prediction 
of a zero cross section near 4.5 GeV, whereas the 
physical A c + D threshold is at 4.15 GeV. Obvi- 
ously the method is inapplicable at low energies, 
which is unfortunate because this regime is of 
greatest interest for RHIC studies of QGP pro- 
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Figure 3: The Kharzeev-Satz J/ip + N total cross 
section Eq.(3.1) and the 1977 SLAC result |J. 

duction. Cross section calculations using this ap- 
proach with more realistic cc wavefunctions are 
currently being carried out by Kopeliovich et al. 




son exchange calculation of J/ip + N — > A c + D 
by Sibirtsev, Tsushima and Thomas [ [l7j finds a 
peak cross section of about 2 mb near y/s = 4.6 
GeV. (An earlier calculation of J/ip + N — > A c + 
D by Haglin assuming D exchange but no 
hadronic form factors, found a somewhat larger 
peak of about 7 mb near y/s — 4.3 GeV.) In com- 
parison, Kharzeev and Satz predict nanobarn- 
scale J/ip + N cross sections 500 MeV above 
threshold, six orders of magnitude smaller! The 
scatter of theoretical predictions for these pro- 
cesses is remarkably wide. 

This discrepancy between theoretical cross 
sections may appear discouraging. One can in- 
stead consider it an opportunity to learn some- 
thing important about QCD: the predictions dif- 
fer because they come from very different as- 
sumptions regarding the scattering mechanism, 
and as they are so far apart, for once we have a 
clear possibility of distinguishing between differ- 
ent hadron-hadron scattering models. 



4. Theory: meson exchange 

Recently, several calculations of charmonium + 
light hadron cross sections have been reported as- 
suming t-channel charmed meson exchange. Of 
course this picture is also problematic, since the 
range of the exchanged charmed meson would be 
only about 1/ATd rs 0.1 fm, and the assumption 
of nonoverlapping hadrons at this separation is 
clearly invalid. (This is the Isgur-Maltman pij 
argument as to why vector meson exchange is un- 
justifed as the source of the short ranged NN core 
interaction.) Nonetheless it is again interesting 
to see what scale of cross section is predicted by 
this type of model, since it might at least incor- 
porate the correct scales and degrees of freedom, 
and it assumes a different scattering mechanism 
from the pQCD color-dipole model advocated by 
Kharzeev and Satz. 

The first such meson exchange calculation, 
due to Matinyan and Miillcr (l^], considered t- 
channel D exchange as the mechanism for the 
reactions J/tp + tt — ► D*D + h.c. and J/ip + p — > 
DD. The cross sections found by this reference 
are shown in Fig. 4 below. Note that 500 MeV 
above threshold these cross sections lie in the 0.5 
to 1 mb range. Similarly, a recent D and D* me- 




Figure 4: The Matinyan-Miiller t-channel meson 
exchange results for J/ip + n and J/ip + p inelastic 
cross sections Jli| . 

Meson exchange calculations of cc + qq dis- 
sociation cross sections have since been reported 
by Lin and Ko Q and Haglin and Gale @|l9|. 
The models differ considerably in detail, due to 
different choices for the diagrams included, the 
effective meson Lagrangian used, and the cou- 
pling constants and form factors assumed. With 
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annihilation- free reactions 25 |2q| . In meson- 
meson scattering there are four distinct diagrams 
(see Fig. 6), each of which has an associated over- 
lap integral of the nonrelativistic quark model 
external meson wavefunctions convolved with the 
interquark Hamiltonian. Constituent interchange 
is forced at Born-order because Hj oc A a • A a 
changes each initial color-singlet qq meson into a 
color octet, which has overlap with color-singlet 
final state mesons only after quark line inter- 
change. The Feynman rules for these diagrams 
were given by Barnes and Swanson |25|. 



°4 4.5 5 5.5 6 6.5 7 

s ,/2 (GeV) 

Figure 5: The f-channel meson exchange cross sec- 
tions for J/ip + N — > A c + D found by Sibirtsev, 
Tsushima and Thomas [jT?|. The smaller contribu- 
tion is from D exchange and the larger is from (non- 
interfering) D* exchange. 

pointlike "hard" vertices, J/tp + it cross sections 
in the 10s of mb not far above threshold are typ- 
ical. After introducing plausible hadronic form 
factors, these are usually reduced to peak val- 
ues of 1-10 mb near y/s w 4.0-4.2 GeV. Unfortu- 
nately there is considerable "guessing" regarding 
hadronic couplings constants, which may be un- 
necessary because these can be calculated using 
well established quark model techniques, for ex- 
ample the 3 Po meson decay model. Similarly, 
hadronic form factors can be derived from the 
3 Po model, presumably with sufficient accuracy 
for these purposes. 

5. Theory: constituent interchange 

We advocate an approach which uses nonrela- 
tivistic quark model wavefunctions and calculates 
these cross sections assuming a constituent in- 
terchange scattering mechanism, driven by the 
Born-order matrix element of the standard quark 
model Hamiltonian. This technique, which has 
no free parameters once quark model wavefunc- 
tions and the interquark Hamiltonian are speci- 
fied, has been shown to compare reasonably well 
with experimental low energy hadron-hadron scat- 
tering data near threshold in a wide range of 




Tl T2 



Figure 6: The four constituent interchange scatter- 
ing diagrams evaluated in the J/ip + qq cross section 
calculation |o[ §1], & The "exchange" is the full 
quark-quark interaction Hamiltonian Hi. 

The first quark model calculation using this 
approach was reported by Martins, Blaschke and 
Quack [|o| , who considered the reactions J/ip + 
7T -> D*D + h.c. and D*D* (The amplitude for 
J/ip + 7r — > DD is zero in the nonrelativistic 
quark model without spin-orbit forces, and has 
been found to be quite weak in a relativized cal- 
culation Martins et al. found that these ex- 
clusive final states have numerically rather simi- 
lar cross sections (except for their different thresh- 
olds), and give a maximum total cross section of 
about 7 mb at yfs ~ 4.1 GeV. A constituent in- 
terchange calculation of J/ip + N and ip' + N 
cross sections using a simplified quark+diquark 
model of the nucleon ]2^] also found several-mb 
peak cross sections not far above threshold. 

Our collaboration recently carried out quite 
similar constituent interchange quark model cal- 
culations (Wong et al. 0]). We used numeri- 
cally determined Coulomb plus linear plus hyper- 



5 



Heavy Quark Physics 5, Dubna, Russia, 6-8 April 2000 



T.Barnes, E.S.Swanson, C.Y.Wong 



fine quark potential model wavefunctions, and 
evaluated the Born-order meson-meson scatter- 
ing amplitude, which is the matrix element of 
the interquark Hamiltonian between scattering 
states with quark interchange (Fig. 6). We in- 
cluded smeared OGE spin-spin, OGE color Cou- 
lomb and linear confining interactions, with pa- 
rameters chosen to give a good fit to the qq and 
cc meson spectra. We find a somewhat smaller 
cross section (peaking at about 1 mb at y/s « 4.0 
GeV) for the sum of these inelastic J/ip + ir re- 
actions. The difference between our work and 
that of Blaschke et al. lies mainly in the treat- 
ment of the confining interaction; for simplicity, 
Blaschke et al. treated confinement as a color- 
independent Gaussian potential that acts only 
between quark and antiquark (hence they include 
only diagrams CI and C2), whereas we used the 
standard A a ■ A a linear confining potential be- 
tween all pairs of constituents. 



■ Tt + J/y->D +DorD+D 

x + J/y -> D + D (a) 

■ Total 71 + J/\\l dissociation cross section 



±1 




Jt + y' ->D + D or D + D 
7i + V|/' ->D + D 

Total 7i + l|f' dissociation cross section 



(b) 



Figure 7: Constituent interchange model predic- 
tions for J/tp + 7r and ip' + n cross sections |2l]] . 
A band of estimated systematic uncertainty is also 
shown. 

We find destructive interference between the 
C and T diagrams, leading to a much reduced 
total cross section. In we also treated ip' + n 
scattering, which involved a simple change to a 
2S cc wavefunction and a change of phase space, 
and found a rather large, ca. 5 mb maximum 



cross section (see Fig. 7). Our cc + qq cross sec- 
tions have their strongest support just a few hun- 
dred MeV in yfs above threshold, since the over- 
lap integrals are damped by the tails of the wave- 
functions at higher energies. 



Total p+7/y dissociation c 




J/\f ->D+D,S m =0 
J/y->D+Dor D+D, S IB =I 
J/\f->D +D , S„=2 
J/\f->D +D',S^=I 
.Aty->o' + o',S,.=» 



Figure 8: Constituent interchange model predic- 
tions for J/tj} + p exothermic reactions as in 
Fig.7. 

We also consider J/tp + p (Fig. 8) and V' + 
p, which are predicted to be much larger near 
threshold for the simple reason that they are 
exothermic; there is a 1/vab divergence in these 
cross sections as we approach threshold. 

Although our scattering amplitudes and cross 
sections are evaluated numerically, it is interest- 
ing that a simple two-parameter function gives 
a useful approximation to our (single channel) 
cross sections. This function is 



a(s) 



max I c 



(5.1) 



where e = yfs — Mq — Mb , and a max (mb) is 
the maximum value of the cross section, at e max 
(MeV) . The threshold exponent p is fixed by the 
angular quantum numbers of the hadrons, and 
is ±1/2 + L^f n (for endothermic/exothermic), 
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where is the lowest angular momentum al- 

lowed for the hnal meson pair CD. As an ex- 
ample, in Fig. 9 we show our numerical results 
for the reaction ry c + tt + — > D + D° and a fit 
using the function (5.1) with p = 1/2, as ap- 
propriate for an S-wave-allowed final state. The 
masses and parameters assumed for this exam- 
ple were M„+ = 0.140 GeV, M Vc = 2.98 GeV, 
M D + = 1.869 GeV, M D o = 1.865 GeV, a s = 0.6, 
b (string tension) = 0.18 GeV 2 , vn u ^ — 0.33 
GeV, m c = 1.6 GeV, and the OGE contact hy- 
perfine smearing distance was 0.25 fm. Those 
are all reasonably well established nonrelativistic 
quark model parameters. The Schrodinger equa- 
tion was solved with these parameters to gen- 
erate numerical wavefunctions, which were then 
used in a nine-dimensional Monte Carlo integra- 
tion to evaluate the scattering amplitudes in the 
CM frame. The Monte Carlo amplitude eval- 
uation in Fig. 9 used AM points for each dia- 
gram at each energy and each final D + (Ct) di- 
rection; amplitudes along three final directions 
were evaluated, which were then projected into 
S-, P- and D-moments to separate the different 
partial waves. Cancellation of diagram sums in 
certain channels as well as evaluation of known 
exact results with SHO wavefunctions provided 
nontrivial checks of the numerical work. 



a[mb] 



sigma(max) = 0.918 mb 
Ecm(max) = 3.802 GeV 



6. Conculsions 

We have reviewed the recent theoretical predic- 
tions and experimental status of the cross sec- 
tions of cc on light hadrons, which is of great 
interest for the interpretation of heavy ion col- 
lisions. There are three scattering mechanisms 
currently being investigated, which are color-di- 
pole pQCD, i-channel meson exchange, and con- 
stituent interchange. The pQCD approach gives 
very small cross sections at low energies, whereas 
the meson exchange and constituent interchange 
models both predict peak cross sections near thres- 
hold of w 1-10 mb. If it is possible to estab- 
lish these cross sections experimentally, we may 
achieve a much better understanding of the mech- 
anisms of low energy hadron-hadron scattering. 
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